Bose-Einstein condensation in tight-binding bands 
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We present a theoretical study of condensation of bosons in tight binding bands corresponding 
to simple cubic, body centered cubic, and face centered cubic lattices. We have analyzed non- 
interacting bosons, weakly interacting bosons using Bogoliubov method, and strongly interacting 
bosons through a renormalized Hamiltonian approach valid for number of bosons per site less than 
or equal to unity. In all the cases studied, we find that bosons in a body centered cubic lattice has 
the highest Bose condensation temperature. The growth of condensate fraction of non-interacting 
bosons is found to be very close to that of free bosons. The interaction partially depletes the 
condensate at zero temperature and close to it, while enhancing it beyond this range below the 
Bose-Einstein condensation temperature. Strong interaction enhances the boson effective mass as 
the band-filling is increased and eventually localizes them to form a Bose-Mott-Hubbard insulator 
for integer filling. 

PACS numbers: 03.75.Lm, 03.75.Nt, 03.75.Hh, 67.40.-w 
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I. INTRODUCTION 

In a many boson system, when the thermal de-Broglie 
wavelength of a particle becomes comparable to the inter- 
particle separation, a condensation in momentum space 
occurs at a finite temperature and a macroscopic num- 
ber of particles occupy the lowest single particle en- 
ergy level and enter in to a phase locked state. This 
phenomenon predicted by Einstein 1 by applying Bose 
statistics 2 to a three dimensional homogeneous system 
of non-interacting atoms in the thermodynamic limit is 
well known as the Bose-Einstein condensation. Though 
it took seventy years to eventually observe Bose-Einstein 
condensation 3 " 6 in metastable, inhomogeneous, finite, 
and three dimensional Bose atom vapors, progress ever 
since has been tremendous. Extensive investigations of 
various aspects of the condensate, a macroscopic quan- 
tum coherent state of atoms, is now a rapidly expanding 
field of research 7-11 . 

One of the interesting recent developments in this field 
is the observation 12 , by Greiner and collaborators, of a 
transition between superfluid and Mott insulating phases 
of bosons in an optical lattice. Experimentalists have 
explored this superfluid-Mott transition and the excita- 
tion spectra in the superfluid 12 ' 13 and Mott insulating 
phases 12 . In these experiments, the condensate is adia- 
batically transferred to a simple cubic optical lattice pro- 
duced by counter propagating laser beams. By changing 
the characteristics of the laser beams, it has been possible 
to achieve great control over t/U, where t is the inter-site 
hopping energy and U the on-site boson-boson interac- 
tion energy. It has been shown that bosons in optical lat- 
tices can be adequately modeled employing a clean Bose- 
Hubbard model 14 . That a band- width controlled tran- 
sition from a superfluid state to a Bose-Mott-Hubbard 
(BMH) insulating state is possible at a critical value of 
t/U for integer number of bosons per site was predicted in 



theoretical studies 14-19 on Bose-Hubbard model. In the 
BMH insulator, the bosons are site localized and the sin- 
gle particle excitation spectrum acquires a gap. Recently, 
superfluid to Mott insulator transition was observed 20 in 
finite one-dimensional optical lattices as well. 



The experimental realization of bosons in optical lat- 
tices provides a microscopic laboratory for the explo- 
ration of collective behavior of quantum many particle 
systems in narrow energy bands with great control on t, 
U, and the number of boson per site (n). There is al- 
ready theoretical studies 21,22 on the possibility of creat- 
ing different types of two dimensional lattices (triangular, 
square, hexagonal, for example). It has been proposed 23 
recently that trimerized optical Kagome lattice can be 
achieved experimentally and that a superfluid-Mott tran- 
sition at fractional filling is possible for bosons in this 
lattice. It is reasonable to expect that three and two di- 
mensional optical lattices of different symmetries will be 
created in near future. Many experimental groups have 
produced three dimensional optical lattices and experi- 
mental studies of Bose condensates in these lattices will 
surely receive increasing attention in the coming years. 
Motivated by such possibilities, we present a theoreti- 
cal study of Bose-Einstein condensation in tight bind- 
ing bands corresponding to simple cubic (sc), body cen- 
tered cubic (bee), and face centered cubic (fee) lattices. 
We have analyzed non-interacting, weakly interacting, 
and strongly interacting bosons. The weakly interacting 
bosons were analyzed using a Bogoliubov type theory 24 . 
For the strongly interacting bosons, we use a renormal- 
ized Hamiltonian valid for n < 1 obtained by projecting 
out on-site multiple occupancies. This analysis is pre- 
sented in the next section and the conclusions are given 
in Sec. III. 
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II. BOSE CONDENSATION IN TIGHT BINDING 
BANDS 

A. Non-interacting bosons 

In this section we discuss the simplest of the three cases 
studied. The Hamiltonian of the non-interacting bosons 
in a tight-binding energy band is: 

H = £[e(k) (1) 

k 

where e(k) is the band-structure corresponding to sc, bee, 
and fee lattices, /i the chemical potential, and is the 
boson creation operator. Confining to nearest neighbor 
Wannier functions overlaps, these band structures when 
lattice constant is set to unity are, 
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£sc(k Xl k yi k z ) — ~2t ^ ' cos(fc^) , 

^L — X 




(3) 

where t is the nearest neighbor boson hopping energy. 
The condensation temperature (Tb) for bosons in these 
bands can be calculated from the boson number equation: 

n = N X N V N Z 51 51 51 e \e(k^k v ,k z )-Alk B T _ I ' ( 4 ) 

y U h h 

r"X " J 'U 

where N s — N x N y N z is the total number of lattice sites, 
ks the Blotzmann constant, and T the temperature. At 
high temperature, the chemical potential is large and neg- 
ative. As the temperature comes down, the chemical po- 
tential raises gradually to eventually hit the bottom of 
the band. Below this temperature there is macroscopic 
occupation of the band bottom and we have a Bose con- 
densate. On further reduction of temperature, the chem- 
ical potential is pinned to the bottom of the band and 
bosons are progressively transfered from excited states in 
to the condensate. All the particles are in the conden- 
sate at absolute zero temperature. Fixing the chemical 
potential at the bottom of the band the solution of the 
number equation gives the Bose condensation tempera- 
ture (Tb). For T > Tb a lower value of [i satisfies the 
number equation, while for T < Tb the RHS of Eq. (4) is 
less than the number of bosons (n), the difference being 
the number of condensate particles (no). We determined 
Tb and no for different lattices as a function of filling and 



temperature. Results of these calculations are shown in 
Figs. 1-5. We find that for bosons in the tight bind- 
ing band corresponding to the bee lattice has the highest 
Bose condensation temperature. Comparing the single 
boson Density Of States (DOS), we find that the band 
structure with smallest DOS near the bottom of the band 
has the highest Tb- The physical reason behind it is that, 
as the temperature is lowered from above the Bose con- 
densation temperature, the bosons are transferred from 
the high energy states to the low energy states following 
the Bose distribution function. To accommodate these 
bosons the chemical potential would touch the bottom 
of the boson band at a higher temperature for a system 
with smaller DOS near the bottom of the band compared 
to a system which has larger DOS there. Consequently 
the Bose condensation temperature for the former system 
would be higher than that for the latter. The growth of 
the condensate fraction for different lattices is shown in 
Figs. (l)-(3), and in Fig (4) we have compared conden- 
sate fraction growth for different lattices. Also shown by 
dotted lines is the condensate fraction for free bosons and 
they are found to be rather close. The variation of Tb 
with n is shown in Fig. (5) which shows an initial fast 
growth and a monotonic increase for higher values of n. 

B. Weakly interacting bosons 

To study Bose condensation of weakly interacting 
Bosons in tight binding bands, we employ the following 
Hamiltonian: 

^ = ^[e(k)- M ]4c k 

k 

+ ^ A rEEE c k+ q c l- q ck ' ck - ( 5 ) 

S k k' q 

Here e(k) is the boson band structure, \x the chemical 
potential, the boson creation operator, U the boson- 
boson repulsive interaction energy taken to be a constant 
for simplicity, and N s the number of lattice sites. In 
this section, we will deal with a range of U such that 
U < 2W where W is the half-band-width. Our aim 
is to get a boson number equation in terms of e(k), U, 
(i, and temperature. One can then calculate the con- 
densate fraction and the transition temperature. To this 
end we follow Bogoliubov theory . In this theory, one 
makes an assumption that the ground state of interact- 
ing bosons is a Bose condensate. Since the lowest sin- 
gle particle state (which is k = in our case of simple 
tight binding bands) has macroscopic occupation (say 
No), we have < ejeo >~< cqc\ >. Then the the op- 
erators c and Co can be treated as complex numbers 
and one gets < c > = < c > = \/No- This com- 
plex number substitution has been recently shown 25 to 
be justified. Clearly we have a two-fluid system consist- 
ing of two subsystems of condensed and non-condensed 
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bosons. There are interactions between particles within 
each subsystem and interaction between particles in the 
two subsystems. In the Bogoliubov approach, to obtain 
second order interaction terms one makes the substitu- 
tion: Cq — > \H% + Cq. In the ground state, the linear n — n 
fluctuation terms terms must vanish and this fixes the 
chemical potential to /i = UriQ + eo, where eo is the min- 
imum of the single particle energy spectrum (which is 
equal to — zt for a bi-partite lattice with co-ordination 
number z), and n = N /N s . Following Bogoliubov ap- 
proach, after a mean-field factorization of the interaction 
term, we obtain for lattices with inversion symmetry the 
following mean-field Hamiltonian: 



H, 



BMF 



-E + J2 



(c k c k + c T _ k c_ k ) 



(6) 



where Eo = —UrioNo/2 and £(k) = e(k) — eo. As men- 
tioned earlier our aim is to get an equation for num- 
ber of particles. This can be obtained from the Green's 
function 26 G(k, u) =« c k ;c k » u) using the relation: 

/oo 
[G(k, uj + in) — G(k, lu — irj)\ f(uj)duj . 
-oo 

(7) 

where /(w) = \/[exp(w /hsT) — 1]. The Heisenberg equa- 
tion of motion for G(k, u) is: 

wG(k,w) = [c k ,c k ]+ << [c k ,i?];4 >> w , (8) 

where H is the Hamiltonian of the system. Using Hbmf, 
we obtain: 

w G(k, u) = 1 + [f (k) + Un }G{k, w) + C/n F(k, w) , (9) 
and the Green's function to which G(k, u) is coupled 

-k' °k 



F(k, w) =<< c^ k ;c k »u> obeys the equation of mo- 
tion: 

uF(k,w) = -[£(k) + [/n ]F(k,w) - J7n G(k,w). (10) 
Solving the above two equations one obtains: 
w + £(k) + Un 



G(k,w) = - 



2£(k) 
1 1 



v w + £(k) w-£7(k) / ' 
where the Bogoliubov quasiparticle energy E k is: 



E(k) = ^ 2 (k) + 2[/n £(k). 



(11) 



(12) 



Note that for the tight binding band dispersions used, 
E(k) is linear in k in the long wave-length limit. Now, 



the number of particles per site (n) is readily obtained 
using Eqs. (7) and (11) to be: 



2N. 



-t 

2N- ^ 

-t 



1 



C(k) + Un 



Ew 



1 



e f3E(k)/k B T _ i 
1 



e -f3E(k)/k B T _ i 



• (13) 



It is useful to look at some limits of the above equation. 
When U — and T — 0, we have no = n which means 
that all the particles are in the condensate at absolute 
zero in the non- interacting limit. When U = and 
T 7^ 0, on obtains: 



n 



1 

n: 



i 



k e k B T 



(14) 
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The second term on the RHS of the above equation is 
the thermal depletion of the condensate. Further, when 
U ^ and T = 0, we get: 



n 



J_j^ / g(k) + Un Q 



E(k) \ 
2E(k) J 



(15) 



in which second term on the RHS is the interaction in- 
duced depletion of the condensate. The interaction has 
twin effects of leading to a modified excitation spectrum 
and to a partial depletion of the condensate. The gap-less 
and linear long-wavelength excitation spectrum is con- 
sistent with experimental measurements 12 on interacting 
bosons in their bose-condensed state in optical lattices. 
In Fig. (6), we have shown the numerical solution of Eq. 
(13) for bee, fee and sc lattices. The condensate fraction 
(no / n) is seen to be gradually suppressed with increasing 
U/W. Note also that at U = 0, all the particles are in 
the condensate. Finally, we consider the case: U ^ and 
T^0. In Fig. (7), we have displayed the variation of the 
condensate fraction as a function of temperature for var- 
ious values of U/W. The interaction partially depletes 
the condensate at zero temperature and close to it, while 
enhancing it beyond this range below the Bose-Einstein 
condensation temperature. We have not plotted these 
curves all the way to Tg since the Bogoliubov approx- 
imation breaks down close to Tg. The variation of Tg 
with n is same as in the case of non-interacting bosons 
as can be seen by setting n = in Eq. (13). At very 
low temperature, thermal depiction is negligible and cor- 
relation induced depletion causes a reduction of no with 
increasing U. At higher temperatures when thermal de- 
pletion is important, U plays another role. The energies 
of the excited states shift to larger values with increasing 
U, consequently the population in the excited states de- 
creases and an enhancement of n occurs with increasing 
U. 

The analysis presented in this section is reasonable pro- 
vided the effect of interaction is perturbative. When the 
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interaction strength increases, there is a possibility for a 
correlation induced localization transition for interacting 
bosons in a narrow band. In the next section, we analyze 
this strongly interacting bosons case. 



C. Strongly interacting bosons 

We first write the Hamiltonian, Eq. (5), in real space. 
Then we have: 



where N and N s are the total number of bosons and 
lattice sites, respectively. In the thermodynamic limit: 



where rii — c\ci. For simplicity, let us confine to the 
case of n < 1. The effect of increasing interaction (U) 
is to make motion of the bosons in the lattice correlated 
so as to avoid multiple occupancy of the sites. In the 
dilute limit, the effect of U is not serious since there 
are enough vacant sites. The effect of U then is promi- 
nent when n is close to unity or to an integer value in 
the general case. In the large U limit, it becomes fa- 
vorable for the bosons to localize on the sites to avoid 
the energy penalty of multiple site occupancy. Quali- 
tatively then, one can see that increasing U increases 
the effective mass (or decreases the band-width) of the 
bosons and eventually drives them, for integer filling, 
to a BMH insulator state. For large U, when the dou- 
ble or multiple occupancy is forbidden one can calculate 
the band-width reduction factor [(/>s(n)] approximately 
following the spirit of Rcnormalizcd Hamiltonian Ap- 
proach (RHA) to fermion Hubbard model 27 ' 28 based on 
Gutzwillcr approximation 29 ' 30 . Within the RHA, the ef- 
fect of projecting out double or multiple occupancies on a 
non-interacting boson wave function is taken into account 
by a a classical renormalization factor which is the ratio 
of the probabilities of the corresponding physical process 
in the projected and unprojected spaces. The probabil- 
ity of a hopping process in the projected space is given 
by n(l — n) for n < 1. This simply implies that the site 
from which hopping takes place must be occupied and 
the target site must be empty in the projected space. In 
the unprojected space, the probability of hopping is just 
equal to the probability of the site from which hopping 
takes place being occupied. The hopping takes place for 
non-interacting bosons irrespective of the target site be- 
ing empty or occupied by any number of bosons. This 
probability may be found out by calculating the num- 
ber of ways a given number of non-interacting bosons is 
distributed in Ns number of lattice sites and (Ns — 1) 
number of lattice sites. The difference would give the 
number of configurations where a particular site is occu- 
pied. Following this route, the probability that a site is 
occupied is obtained as: 



p(N;N s ) = l 



(N + N S -2)\(N S -1)! 

(n s - 2)\(n + n s - iy: 



p(n) 



l + n 



(18) 



So, the hopping probability is just p(n). The above equa- 
tion is valid for any n. Now, in the strongly correlated 
state (large U limit), confining to the case of n < 1, the 
hopping probability is n (1 — n). Hence the 4>b (n) is ob- 
tained to be, 



(j> B {n) = l- n 2 . 



(19) 



The above equation is valid only for n < 1 and in the 
large U limit. It may be interesting for the reader to 
note that, in the Fermion case, 0i?(n) = 2(1 — n)/(2 — n) 
which has been used in the studies 31 of superconductiv- 
ity in strong coupling fermion Hubbard model of high 
temperature superconductors. Since in the large U limit, 
double or higher site-occupancies are forbidden, one can 
write a renormalizcd Hamiltonian, valid for n < 1, for 
strongly correlated bosons as : 



H sc = ^2[<pB(n)e(k) - ^]c\.c^. 



(20) 



(17) 



The above H 8C is clearly the Hamiltonian of non- 
interacting bosons in a narrow band which has under- 
gone a strong correlation induced filling-dependent band- 
narrowing. One can see that as n increases from to 1, 
the boson effective mass increases to eventually diverge 
at n = 1 and a BMH insulator obtains. We do admit 
that there are limitations to this renormalized Hamilto- 
nian approach. While it has a merit that detailed band 
structure information can be incorporated in the Bose 
condensation temperature calculation, it has the demerit 
that we have to restrict ourselves to large U and n < 1. 
The above Hamiltonian is valid for U/W > (U/W) c 
where (U/W) c is the critical value for transition into the 
Mott insulating phase for n = 1. Fixing a precise lower 
limit on (U /W) is not possible in the absence of either 
exact analytical or numerical solution of three dimen- 
sional Bose-Hubbard model. It may be mentioned that 
in this large U limit and for n < 1, the Bose-Hubbard 
model is reduced to the lattice Tonks (hard-core boson) 
gas which is different from the classical gas of elastic hard 
spheres investigated by Tonks 32 . Our results imply then 
that the effective mass of bosons in a lattice Tonks gas in 
a narrow energy band is strongly band-filling dependent. 
The variation of bose condensation temperatures with n 
for bosons with correlation-induced renormalized energy 
bands corresponding to sc, bee, and fee lattices are dis- 
played in Fig. (8). In the region between each curve and 
the n-axis, the bosons are in their bose-condensed state, 
and above each curve they are in their normal state (ex- 
cept at n = 1 and below fc^T m U). The variation of 
Tb is a resultant of the combined effects of increasing 
density and increasing effective mass since Tb is propor- 
tional to n/m*. Beyond around twenty percent filling, 
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the increasing effective mass over-compensates the effect 
of increasing n and pulis down the growth of Tb eventu- 
ally driving it to zero at n = 1 at which density one has 
a BMH insulator. 

It is of some interest to make a comparison be- 
tween Mott-Hubbard (MH) metal-insulator transitions 
observed in fermion systems in condensed matter physics. 
For a half- filled narrow band of fermions, one way to in- 
duce a MH insulator to metal transition is by reducing 
the ratio of Coulomb repulsion (U) to the band-width 
(2W). This was experimentally achieved 33 in V2O3 by 
application of pressure. This then is a band-width con- 
trolled MH transition 34 . Another way to induce a MH 
insulator to metal transition is to start with a Mott in- 
sulator and reduce the band-filling which would then be 
filling-controlled MH transition. This was achieved 35 in 
La x Sri- x TiO?,. For a simple Gutzwiller approximation 
based analysis of the properties of this material see Ref. 
36. Now, the Mott transition observed 12 in boson sys- 
tems in optical lattices is the band-width controlled one. 
It would be interesting to look for filling-controlled Mott 
transition in boson systems in optical lattices which may 
be possible by starting with the BMH insulator and flip- 
ping a few atoms out of the trap. 



III. CONCLUSIONS 



In this paper we presented a theoretical study of 
condensation of bosons in tight binding bands corre- 
sponding to sc, bec, and fee lattices. We analyzed 
condensation temperature and condensate fraction of 
non-interacting bosons, weakly interacting bosons us- 
ing Bogoliubov method, and strongly interacting bosons 
through a renormalized Hamiltonian approach (limited 
to n < 1) capable of incorporating the detailed boson 
band structures. In all the cases studied, we find that 
bosons in a tight binding band corresponding to a bec lat- 
tice has the highest Bose condensation temperature. The 
growth of condensate fraction of non-interacting bosons 
is found to be very close to that of free bosons. In the 
case of weakly interacting bosons, the interaction par- 
tially depletes the condensate at zero temperature and 
close to it, while enhancing it beyond this range below 
the Bose-Einstein condensation temperature. Strong in- 
teraction enhances the boson effective mass as the band- 
filling is increased and eventually localizes them to form a 
Bose- Mott-Hubbard insulator for n = 1. In the strongly 
interacting bosons case, we found that all bosons are in 
the condensate at absolute zero temperature. We also 
pointed out a possibility of a filling-controlled BMH tran- 
sition for bosons in optical lattices. 
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FIG. 2. The same as in Fig. 1 for bosons in a bee lattice 
for n = 0.8 (top), n = 0.6 (middle), n = 0.4 (bottom). 
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FIG. 1. The varation of condensate fraction with temper- 
ature for bosons in sc lattice for n — 0.8 (top), n = 0.6 
(middle), n = 0.4 (bottom). In this and other figures W is 
the half-band-width. 
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FIG. 3. The same as in Fig. 1 for bosons in a fee lattice 
for n = 0.8 (top), n = 0.6 (middle), n = 0.4 (bottom). 
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FIG. 4. The variation of condensate fraction (for n = 0.6) 
for bcc (top), fee (middle), and sc (bottom) lattices. The dots 
are plots of 1 - (T/T c )'^ 2 . 



FIG. 6. The variation of condensate fraction (for T — 
and n = 0.4) with U/W for bcc (top), fee (middle), and sc 
(bottom) lattices. 
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FIG. 7. Condensate fraction (for n = 0.4) vs temperature 
for various values of U/W (shown on the curves) and for sc 




bcc (top), fee (middle), and sc (bottom) lattices. 
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FIG. 8. Bose condensation temperature (of strongly inter- 
acting bosons) vs n for bcc (top), fee (middle), and sc (bot- 
tom) lattices. The dots are for non-interacting bosons in a 
bcc lattice. 



